Mon. Not. R. Astron. Soc. 000, 000-000 (0000) 



Printed 1 February 2008 



(MN LMftX style file v2.2) 



■ 

o ■ 
o : 

(N . 

Oh: 



A POSSIBLE STELLAR METALLIC ENHANCEMENT IN 
POST-T TAURI STARS BY A PLANETESIMAL 
BOMBARDMENT 

O. C. Winter 1 *, R. de la Reza 2 , R. C. Domingos 1 , L. A. G. Boldrin 1 and C. Chavero 5 

1 Sao Paulo State University - UNESP, Grupo de Din&mica Orbital & Planetologia, Guaratinguetd, CP 205, CEP 12.516-410, Brazil. 

2 Observatdrio Nacional, Rio de Janeiro, RJ, Brazil 

1 February 2008 

ABSTRACT 

The photospheres of stars hosting planets have larger metallicity than stars lacking 
planets. This could be the result of a metallic star contamination produced by the 
bombarding of hydrogen deficient solid bodies. In the present work we study the 
possibility of an earlier metal enrichment of the photospheres by means of impacting 
planetesimals during the first 20-30Myr. Here we explore this contamination process by 
simulating the interactions of an inward migrating planet with a disc of planetesimal 
interior to its orbit. The results show the percentage of planetesimals that fall on the 
star. We identified the dependence of the planet's eccentricity (e p ) and time scale of 
migration (r) on the rate of infalling planetesimals. For very fast migrations (r = 10 2 yr 
and t = 10 3 yr) there is no capture in mean motion resonances, independently of 
the value of e p . Then, due to the planet's migration the planetesimals suffer close 
approaches with the planet and more than 80% of them are ejected from the system. 
For slow migrations (r = 10 5 yr and r = 10 6 yr) the percentage of collisions with the 
planet decrease with the increase of the planet's eccentricity. For e p — and e p = 0.1 
most of the planetesimals were captured in the 2:1 resonance and more than 65% 
of them collided with the star. Whereas migration of a Jupiter mass planet to very 
short pericentric distances requires unrealistic high disc masses, these requirements are 
much smaller for smaller migrating planets. Our simulations for a slowly migrating 0.1 
Mjupitcr planet, even demanding a possible primitive disc three times more massive 
than a primitive solar nebula, produces maximum [Fe/H] enrichments of the order 
of 0.18 dex. These calculations open possibilities to explain hot Jupiters exoplanets 
metallicities. 
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1 INTRODUCTION 

A particular interesting problem concerning exoplanetary 
science is metallicity. Even if the near 200 extrasolar giant 
planets that have been detected by transit and mainly ra- 
dial velocity methods are concentrated around in general 
nearby solar type stars, a clear metallicity correlation has 
been established. In fact, stars with planets (SWP) appear 
to be more metal rich than stars without planets. A general 
metallicity shift of 0.2 dex in [Fe/H] is characteristic for 
SWP (Santos et al. 2005). In particular it is interesting to 
note that stars with hot Jupiters (with orbital periods less 
than about 5 days) appear to be more metal rich than the 
mean of all SWP (Sozzetti et al. 2004; Fischer & Valenti 



* ocwinter@feg.unesp.br 



2005; Butler et al. 2006). These results are concentrated in 
surveys and analyses of non-metal poor stars and for spec- 
tral types F (young stars), G and K (old stars). Two impor- 
tant review articles devoted to the chemical composition of 
SWP, indicating also properties of other elements than Iron, 
appeared with results obtained up to 2003 (Gonzalez 2003) 
and for the period 2003-2006 in Gonzalez (2006). 

What is the physical mechanism behind these metallic- 
ity relations ? Two main different interpretations have been 
proposed in the literature, nevertheless, any of them has 
produced a satisfying explanation. On the one hand, is a 
primordial condition in which metal rich SWP would have 
been formed in metal rich clouds. On the other hand, there 
is the external accretion mechanism. Here the stellar sur- 
faces, containing poor or normal metal abundances enhance 
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their metallicities by the injection of solid metal rich mate- 
rial depleted in H and He. 

These two mechanisms will form, metallically speaking, 
two different types of stars. The initial or primordial mecha- 
nism, produce entire metal rich stars, from their centres up 
to their atmospherical layers. This is not the case for the 
accretion scenario, which can modify the metal content of 
the external layers only, at least for those stars with shallow 
convective layers. The primordial mechanism is only con- 
sidered in the literature when any argument in favor of the 
self-enrichment is discarded. Anyway, the primordial mech- 
anism necessitates a highly inhomogeneous metallicity dis- 
tribution of the interstellar matter, not detected at present. 
A powerful argument in favor of the primordial case, would 
consist in finding that the centres of SWP by asteroseismo- 
logical methods, are also metal rich. The only star studied 
with this technique was /i Arae (Bazot & Vauclair, 2004; 
Bazot et al. 2005). 

Independent of this, the accretion or self-enrichment 
mechanism can continue to be explored and this is the pur- 
pose of the present work. We study here the possibilities of 
bombardment of stars atmospheres by their discs planetesi- 
mals during a certain phase of the pre-main sequence (PMS) 
evolution. We hope that this will provide some clues respect 
to metallic enrichment. 

In the next section we discuss some aspects of the evolu- 
tion of PMS discs. The numerical model is presented in Sec- 
tion 3. The results from our simulations are given in section 
4. In Section 5 some problems related to the mass of discs are 
discussed. Section 6 is devoted to the metal enhancements 
and finally in Section 7 we present our conclusions. 



2 PRE-MAIN SEQUENCE DISC EVOLUTION 

Several stages, some of them badly known, are involved in 
discs evolution in the PMS. If discs during the T Tauri phase 
(before 5 Myr) are formed by a relative homogeneous distri- 
bution of gas and dust, during the next post- T Tauri phase 
the situation is very different. In fact, in post- T Tauri stars 
(PTTS) disc gas and dust components follow quite different 
histories. Gas is lost in a canonical period of 10 Myr (see 
for instance Sicilia-Aguilar et al. 2005). However. PTTS as- 
sociations belonging to the older subgroups of the Sco-Cen 
OB association appear to retain their disc gas up to ages of 
16 Myr (Pinzon et al. 2007). In any case, during the first 
30 Myr, the fine dust content is gradually agglutinated into 
kilometric sized rocky or iced bodies (planetesimals). Even- 
tually, in the outer parts of the discs, conditions are fulfilled 
for the formation of a core of 10 - 20 MEarth capable to form 
a giant planet, attracting the gas when this is still available. 

When gas is almost absent, dust evolution can be more 
complex because energetic collisions between planetesimals, 
produce a new generation fine dust forming what is called 
a Debris Disc (DD) which can be maintained up to very 
large ages, even Gyr. At present, our knowledge of planetary 
formation in DD is confused. 

Considering the above discussion, it is possible to con- 
ceive the existence of discs with ages between 10 and 30 
Myr in which an external planet coexist with a sea of in- 
ternal planetesimals. If the total mass of planetesimals is 
of the same order of the planet, interactions between these 



two components will provoke an internal migration of the 
planet. This because the planet is continously loosing en- 
ergy and the excess energy is being used to disperse the 
planetesimals (Murray et al. 1998). Since their work, plane- 
tary migration has become an important part in the study of 
planetary systems. Some recent lectures or reviews on this, 
are available in the literature. For extra solar systems in Ar- 
timowicz (2006) and for the solar system in Levison et al. 
(2007). More particularly, several recent studies appeared 
using N-body simulations in order to study the conditions 
of survival (or formation) of terrestrial type planets due to 
migration of a giant planet (Lufkin, Richardson & Mundy 
2006; Fogg & Nelson 2005; Raymond et al. 2006; Armitage 
2003; Mandell & Sigurdsson 2003). In respect to N-body sim- 
ulations of planetesimals and a migrating planet directed to 
the metallicity enhancement problem, the only work to our 
knowledge is that of Quillen & Holman (2000) (hereafter 
QH00). In principle our work here is an extension of QH00. 



3 NUMERICAL MODEL 

In this work, we numerically integrated planar systems com- 
posed by a star, a planet and one thousand planetesimals. 
It was assumed that the star has the same mass as our Sun 
and the planet the same mass as Jupiter. The planet is ini- 
tially with semi-major axis a p = 5AU and eccentricity e p . 
The planet is forced to migrate inward up to a v = 0.01AU, 
with constant speed in a timescale r. The planetesimals are 
considered to be massless particles that are initially dis- 
tributed on random circular orbits with semi-major axis 
1 < a < 4AU. It is important to note that this migration 
is not self-consistent in that the bodies ejected do not cause 
migration. 

The numerical simulations were performed using a code 
from SWIFT (Levison & Duncan, 1994). The integrator is 
based on the MVS method developed by the Wisdom & 
Holman (1991). The output of the simulations give the tem- 
poral evolution of the semi-major axis and eccentricity of 
the planetesimals. The integration for each particle stopped 
when one of the three following conditions happened: 

a - collision with the star: when the planetesimal gets 
closer than 0.01AU from the star; 

b - collision with the planet: when the planetesimal gets 
closer than two Jupiter's radius from the planet; 

c - ejection from the system: when the planetesimal 
reaches more than 50AU from the star. 

The value of the time scale, r, is defined by the process 
responsible for the planet's migration Two main process are 
usually discussed in the literature. One is the migration of 
the planet due to the gravitational interaction with a disc 
of planetesimals (Murray et al. 1998). The other is the mi- 
gration of a planet embedded in a gaseous disc (Artymow- 
icz 2005). In the first case the time scale is of the order 
of 10 5 — 10 7 yr, while in the second case it is much faster, 
10 2 - 10 3 yr. In our studies we tried to cover the whole 
spectrum of time scales. 
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4 NUMERICAL SIMULATIONS 

In order to study the influence of the planet's eccentricity on 
the spreading of the planetesimals we performed integrations 
for 0.0 ^ e p ^ 0.5 with Ae p = 0.1. In Figures 1 to 6 are pre- 
sented representative snapshots of the temporal evolution 
of the semi-major axis and eccentricity of the planetesimals, 
where r = 10 . There is a strong connection between the 
orbital evolution of the planetesimals and mean motion res- 
onances between the planetesimals and the planet. In each 
plot are indicated the location of the semi-major axis for the 
main mean motion resonances. 

In Figure 1 are presented axe diagrams for the sim- 
ulation with e p = 0. The diagrams show the state of the 
planetesimals at t=0, at t=l x 10 3 yr, at t=2 x 10 4 yr and 
at t=1.5 x 10 'yr. All planetesimals start with circular orbits 
and at some stage are captured in a resonance. Then, their 
eccentricities start to growth until they reach a very high 
value and be removed from the system by collision or ejec- 
tion. In this case (e p = 0) the planetesimals are captured 
in one of the following resonances: 3:2, 5:3 and 2:1. The 2:1 
resonances is the one that plays the main role in the orbital 
evolution of the planetesimals. 

In Figures 2 to 4 are presented axe diagrams for the 
simulations with e p = 0.1, 0.3, and 0.5, respectively. As the 
value of the planet's eccentricity increases more mean mo- 
tion resonances, located further from the planet, become im- 
portant for the dynamics of the planetesimals. While those 
closer to the planet become less important. In the case of 
e p = 0.1 the resonances 5:2 and 3:1 are also important, but 
the resonance 2:1 is still the one that plays the main role in 
the orbital evolution of the planetesimals. For e p = 0.2 the 
resonance 4:1 is also important while the resonances 3:2 and 
5:3 contribute very little in the orbital evolution of the plan- 
etesimals. For e p — 0.3, 0.4 and 0.5 other resonances become 
also important (7:2, 5:1, 6:1, ...). Then the orbital evolution 
of the disc of planetesimals is affected by several resonances 
and is not dominated by a particular one. In these cases the 
planetesimals of the outer part of the disc (a>0.65AU) are 
spread without being captured in resonance. 

Along each numerical simulation the number of plan- 
etesimals goes decreasing with time. The temporal evolution 
of the remaining percentage of planetesimals in the system 
is shown in Figure 5. Each plot of Figure 5 corresponds to 
the simulation with a given e p , where m p = Mj up iter and 
r = 10 6 yr. The plots show that with the increase of the 
planet's eccentricity the number of remaining planetesimals 
in the system decreases faster. 

In order to study the influence of the migration speed on 
the spreading of the planetesimals we performed integrations 
for different values of timescales, from r = 10 2 yr to r = 
10 6 yr with At = 10 yr. Simulations for r = 10 7 yr are 
computationally too expensive. We made a few simulations 
with r = 10 7 yr and the results are not much different from 
those for r = 10 6 yr. 

In Figure 6 we present the statistics of the results in 
terms of percentage of collisions with the star, with the 
planet and ejections from the system. For slow migrations 
(r = 10 5 yr and r = 10 6 yr) the percentage of collisions 
with the planet decrease with the increase of the planet's 
eccentricity. For e p = and e p = 0.1 most of the planetes- 
imals were captured in the 2:1 resonance and more than 



65% of them collided with the star. For very fast migrations 
(r = 10 2 yr and r = 10 3 yr) there is no capture in mean 
motion resonances, independently of the value of e p . Then, 
due to the planet's migration the planetesimals suffer close 
approaches with the planet and more than 80% of them are 
ejected from the system. 

Quillen (2006) developed analytical predictions of cap- 
ture in resonances. His results show that raising e p above 
a certain limit, en m , may prevent the resonance capture. 
He also predicts that the capture probability is reduced for 
faster migration. Our results corroborate his predictions. 
The capture probability for resonance also depends on the 
planets mass (Quillen, 2006), which will be changed in the 
next section. 



5 MASS PROBLEMS 

The minimum mass in scattering planetesimals necessary to 
produce the migration of a planet can be given by (Adams 

6 Laughlin 2003) 

M — m p ln(oo/a/) 

where m p is the mass of the planet, ao the initial semi-major 
axis of the planet and a/ its final semi-major axis. From our 
simulations , ao = 5 AU, a/ = 0.01 AU and m p equal to 
one Jupiter we obtain that M is equal to about 6 Mj up itcr. 
As we will see after the migration of a Jupiter mass planet 
up to the pericentre distance of 0.01 AU (2.15 solar radius) 
demands extremely large disc masses. However, this is not 
the case for a 0.1 Mj up j t er migrating planet. Davis (2005) 
proposed a cumulative mass growth in the primitive solar 
nebula model such that the total mass of the disc up to 
Jupiters position would be about 3838 MEarth- Adopting a 
canonical, gas to dust ratio of 100, we expect only 1% in 
the mass of planetesimals. Then the mass of planetesimals 
in our adopted disc model from Davis (2005) would be 38 
AfEarth- Our simulations show for a planet with Jupiters 
mass, in near circular orbit, migration from 5 AU to 0.01 
AU in a time scale of r > 10 5 yr, that about 80% of the 
planetesimals collide with the star. However, the minimum 
mass required for a Jupiter mass planet to migrate up to 
that pericentre distance is, as we see above, of the order of 6 
Afjupiter. This represent for ( Mj up jter= 318 MEarth) a disc 
mass of 1908 MEarth of planetary mass, which represent 50 
times our adopted primitive solar nebula. This is unrealistic. 
The situation is much better for a migrating planet of 0.1 
Mjupiter with a required mass of 190 MEarth- 

As an example, lets consider which will be the re- 
quirements for the maximum concentration of observed hot 
jupiters exoplanets at 0.05 AU. In this case, The necessary 
minimum mass is 1464 MEarth which is a factor 38 larger 
than our primitive solar adopted mass and represent a disc 
of 0.4 solar mass, also unrealistic. For the case of 0.1 M.j up itor 
planet, for the same migration distance requires a disc (up 
to Jupiter distance) of 0.04 solar mass which could be ap- 
propriate. 

Considering this, we search for possibilities for a slow 
migration ( r = 10° - 10 yr) of a 0.1 Mj up iter mass planet. 
The simulations results for this case are shown in Table 2. 
Here, the planet is set to migrate from 5 AU up to 0.01 AU 
in principle with r= 10 6 yr. We performed simulations for 
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Figure 1. : Temporal evolution of the semi-major axis and eccentricity of the planetesimals for the simulation with m p = Mj up i ter , 
r = 10 6 yr and e v = 0. The axe diagrams show the state of the planetesimals at t=0, at t=l x 10 3 yr, at t=2 x 10 4 yr and at t=1.5 x 10 5 yr. 
In each plot are indicated the location of the semi-major axis for the main mean motion resonances. 
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Figure 2. : Temporal evolution of the semi-major axis and eccentricity of the planetesimals for the simulation with m v = Mj up j t er, 
r = 10 6 yr and e p = 0.1. The axe diagrams show the state of the planetesimals at t=l X 10 3 yr, at t=2 X 10 4 yr, at t=8 X 10 4 yr and at 
t—1.2 x 10 5 yr. In each plot are indicated the location of the semi-major axis for the main mean motion resonances. 
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Figure 3. : Temporal evolution of the semi-major axis and eccentricity of the planetesimals for the simulation with m p = Mj up i ter , 
r = 10 6 yr and e p = 0.3. The axe diagrams show the state of the planetesimals at t=3 x 10 3 yr, at t—2.2 x 10 4 yr, at t=3.8 x 10 4 yr and 
at t=6.2 x 10 4 yr. In each plot are indicated the location of the semi-major axis for the main mean motion resonances. 
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Figure 4. : Temporal evolution of the semi-major axis and eccentricity of the planetesimals for the simulation with m v = Mj up it e r, 
r = 10 6 yr and e v = 0.5. The axe diagrams show the state of the planetesimals at t=l X 10 3 yr, at t=6 X 10 3 yr, at t=1.6 X 10 4 yr and 
at t—2A x 10 4 yr. In each plot are indicated the location of the semi-major axis for the main mean motion resonances. 
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Figure 5. : Temporal evolution of the remaining planetesimals in the system. These are the results from the simulations also presented 
in Figures 1 to 4, where m v = Mj up it e r and r = 10 6 yr . The eccentricity of the planet is: a) e p = 0, b) e p = 0.1, c) e v = 0.2, d) e p = 0.3, 
e) e p = 0.4 and f) e p = 0.5. 



Table 1. Summary of the results for the simulations with m p = Mj up it er and e p = 0, its is given the percentage of 
collisions of planetesimals with the star as a function of the planet's migration time scale. It is also shown comments 
on the problems related with each case. 



r(yr) Migration Model Star Collision (%) Problem 

10 2 Gas >5 low star collision rate 

10 3 Gas >5 low star collision rate 

10 4 - 20 — 

10 5 Planetesimals 70 disc too massive 

10 6 Planetesimals 70 disc too massive 



e p = to e p = 0.5 with Ae p — 0.1. In columns 5 and 6 are 
given the time and the value of the semi-major axis of the 
planet when the last planetesimal is removed from the sys- 
tem. For e p = 0.1 (Figure 7) the last planetesimal is removed 
at t = 1.71 x 10 5 yr, when the planet has a f = 1.72 AU. In 



order to produce such migration the mass of ejected plan- 
etesimals (Equation 1) has to be of the order of 34 MEarth. 
We found that in this simulation 50% of the planetesimals 
collide with the star and about 45% are ejected (Table 2). 
Therefore, in such case a disc of planetesimals with mass of 



STELLAR METALLIC ENHANCEMENT BY A PLANETESLMAL BOMBARDMENT 7 



100 



03 

£= 

J5 

CL 




100 



t = 10" yrs 




t = 10 yrs 




t = 10 yrs 




100 



t = 10 6 yrs 




collision with the star 



-•— collision with the planet 



escape 



Figure 6. : Planetesimals end state for the simulations with m p = Mj up itor- In each plot is given the percentage of planetesimals that 
collided with the star (triangle), that collided with the planet (circle) and that escaped (square). The values are given as a function of 
the planet's eccentricity. On the top of each plot is shown the corresponding time scale for the migration of the planet. 
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about 70 AfEarth would not be too massive for the present 
models (Davis 2005) and would produce the planets migra- 
tion required to generate the metallicity expected. 



6 METAL ENHANCEMENT IN PTTS 

Any model using injection of rocky planetesimals as a way to 
enhance stellar metallicities of young stars, requires a com- 
promise between the epochs of existence of planetesimals, 
the time of planet migration (at least for a case of a single 
planet as is the case here) and the epoch of diminishing of 
the surface convective layer. This last condition is specially 
important in order that the metal enhancement be efficient. 
In fact, very young stars as T Tauri, are characterized by 
their large convection zones and these must be stabilized to 
their minimum configuration. In that case, the new injected 
material is trapped and not diluted as is the case in stars 
with large convective zones. 

Due to a very rapid decrease of convective envelopes 
of solar type stars, this mentioned compromise appears be- 
tween 20 and 30 Myr. This "window of opportunity" as is 
clearly shown in Figure 5 in Ford, Rasio & Sills (1999) rep- 
resent the right time to enhance the metal content, that in 
principle is maintained and observed today. A possible sig- 
nature of this metallic enrichment has been detected in de 
la Reza et al. (2004). Here a collection of G and K stars 
belonging to two coeval PTTS associations with ages of 20 
and 30 Myr showed G stars with larger Fe abundances in 
respect to the lower mass (and larger convection zones) K 
type stars. In our scenario, as is also the case of QH00, the 
metal enrichment can be produced by particles directly im- 
pinging in the stars or by the grazing planetesimals with 
very large eccentricities. In our simulation we obtain , for 
the low migrating rates and depending on the approaching 
distances to the star, 10 3 up to 10 5 close passages. 

QH00 considered that if rocky bodies are larger than 
one kilometer, they will survive vaporization, other possible 
destiny of particles is disruption of strengthless material. 
Sridhar & Tremaine (1992) (see also Asphaug & Benz 1996) 
found that inviscid planetesimals in parabolic orbits will be 
disrupted by the planet, or in our case by the star, if their 
pericentric distances are less than 1.69 R s {p s / pp) 1 ^ 3 • Here, 
R s is the stellar radius, p s the stellar density and p p the den- 
sity of the planetesimal. After disruption, the planetesimals 
became a kind of "dust" maintaining their same densities. 

In view of this and as an example, we can explore which 
will be the disruption distances for three types of represen- 
tative stars. An A type star with (two solar mass, two solar 
radius), a solar G star, and a K star (half solar mass and 
radius). Their respective densities being 0.35, 141 and 5.6 
g/cc. The corresponding disruption distances will be 0.8, 
1.31 and 2.1 stellar radius. Taking into account that in our 
simulations, we consider that an impact is defined when a 
particle is closer than 0.01 AU (2.15 solar radius), we have 
that for our A star with 2 solar radius, practically all par- 
ticles penetrate the star and "disruption" will occur inside 
the star (of course this is not valid because the physical sit- 
uation inside the star must be different). For a solar type 
star, particles with a pericentric distances less than 1.3 so- 
lar radius will be disrupted. For K stars particles already at 
distances less than 2 stellar radius will be disrupted. In any 



case, it is expected that these particles will contribute to the 
metal enrichment of the star. 

Now, lets consider which will be the possible metal en- 
hancements for the most realistic case in which a planet of 
0.1 Mjupitor is migrating. In Table 2 are displayed, for dif- 
ferent planets eccentricities, the minimum masses (Col. 7) 
following relation 1, required to terminate the migration at 
the pericentric distances indicated in (Col. 6). The impact- 
ing masses are shown in (Col 8). In case of high penetration 
efficiency of particles, we obtain that the best cases are rep- 
resented by planet eccentricities 0.0, 0.2 and 0.3. Adopting 
the Fe enrichment calculated curves in Figure 5 by Ford et 
al.(1999), we obtain for the impacting masses, [Fe/H] abun- 
dances at 20 and 30 Myr, of the order of 0.18 dex for e p = 
0.0 and e p = 0.2 and about 0.15 dex for e p = 0.3. Because for 
our simulations shown in Table 2, we need planetary max- 
imum masses of the order of 100 M^rth- Reminding that 
our primitive solar disc contains about 38 ME ar th of solid 
matter, we conclude that primitive discs to produce these 
enrichments must be larger than a factor 3 than the solar 
one. Larger factors, nevertheless possible to exist, are neces- 
sary to explain the observed hot Jupiters exoplanets at 0.05 
AU, with masses below one Jupiter mass, even with sightly 
larger Fe abundances. 



7 MAINTENANCE OF THE METALLIC 
ENHANCEMENT 

The problem of the long term preservation of the metallic 
excess in the surface stellar layers obtained by accretion is 
a difficult one. Let us first consider the situation concern- 
ing the convective layers. As mentioned in this work, the 
accretion of solid matter devoid of hydrogen refers only to 
the contamination of the surface convective layers. In this 
way, the sizes of these layers are of a prime importance. 
The present majority of SWP belongs to FGK types stars 
and future results on large searches on A and M type stars 
will bring more light on this problem. Due to the impor- 
tance of M dwarfs as being an important part of the stellar 
content of the Galaxy, planet searches are active for these 
stars. Nevertheless, only few unities of M dwarfs planet hosts 
have been discovered, for example, GJ876 (a multiple planet 
host), GJ 581 and GJ 436 (only with candidate planetary 
companions), even few, they are producing interesting and 
instructive results. Recent detailed spectroscopic analyses of 
these stars by Bean, Benedict & Endl (2006) and another 
using photometry (Bonfils et al. 2005), indicates that all 
these dwarfs are metal deficient. At present it is not clear 
if these results are a consequence of a special metallic dis- 
tribution of the M stars candidates for the planet searches, 
but independently of this, we can consider here that due to 
their very large diluting convective zones, no supersolar M 
dwarfs planets hosts are expected if the accretion mecha- 
nism is working. In this way, an eventual future discovery 
of a supersolar M dwarf planet host will contribute to the 
validity of the primordial mechanism as mentioned in the 
introduction section. 

Concerning the question if any metallic enhancement 
produced in the first 20 - 30 Myr will be maintained up to 
Solar ages or larger, we do not have a clear answer yet. In 
fact, even if the convective zones are relatively constant from 



STELLAR METALLLC ENHANCEMENT BY A PLANETESLMAL BOMBARDMENT 9 



t = 1 .0 x 1 0° yrs 




1 

0.8 
0.6 
0.4 
0.2 




t = 6.0x 10 J yrs 



6 1 4? 2:1 3:2 

■ 5:1}7:2< 





Figure 7. : Temporal evolution of the semi-major axis and eccentricity of the planetesimals for the simulation with m p = 0.1Mj up i tcr , 
r = 10 6 yr and e v = 0.1. The axe diagrams show the state of the planetesimals at t=l x 10 3 yr, at t=6 x 10 3 yr, at t=3.6 x 10 4 yr and 
at t=6 x 10 4 yr. In each plot are indicated the location of the semi-major axis for the main mean motion resonances. 



Table 2. Simulations for a small planet (m p = 0.1Mj up i tor ) migrating from 5AU to 0.01AU in a time scale r = 10 f, yr. In 
columns 5 and 6 are given the time and the value of the semi-major axis of the planet when the last planetesimal is 
removed from the system. In column 7 is given the minimum masses (Eq. (1)) required to terminate the migration at 
the pericentric distances indicated in column 6. The impacting masses are shown in column 8. 



e p Star Collision (%) Planet Collision (%) Ejected (%) <iast(xl0 5 )yr a p _i a st(AU) Af(M Earth ) Impact. Mass (M Earth ) 



0.0 68.5 3.1 28.4 2.52 1.04 50 34 

0.1 50.0 5.5 44.5 1.71 1.72 34 17 

0.2 43.6 7.3 49.1 3.70 0.50 73 32 

0.3 35.7 5.8 58.5 5.29 0.19 104 37 

0.4 25.8 6.3 67.9 1.85 1.59 36 9 

0.5 20.3 6.1 73.6 1.65 1.79 33 7 



10 8 up to 10 9 yr (even diminishing somewhat near 10 10 yr 
as can be seen in Fig. 5 of Ford et al. 1999) other mech- 
anisms as long term diffusion can be acting. As suggested 
by Gonzales (2006), one interesting way to investigate the 
maintenance of the excess metallicity created by accretion, 
consists in doing planet searches on specific open clusters. 
For example, in M67 with an age and metallicity similar to 
that of the Sun. In these kind of researches with open clus- 
ters, where ages and metallicities are no more variables, any 
SWP which would be enriched during the pre-main sequence 
would present smaller metallicities during the sub-giant and 
giant phases, testing this way the accretion mechanism. 

At present, the only cluster surveyed for planet searches 
is the Hyades, due to its favorable conditions as proximity 
and high metal abundance. Even if no planets have been de- 
tected among their dwarf stars (Cochran, Hatzes & Paulson 



2002, Paulson, Cochran & Hatzes 2004), the first planet (a 
massive one), ever discovered in a cluster, has been recently 
detected around a genuine member of the Hyades (Sato et 
al. 2007). The host star is the classified K0III clump gi- 
ant Eps Tau (also massive) with a measured metallicity of 
[Fe/H] = 0.17 ± 0.04, a value similar to that of the Hyades 
cluster ([Fe/H] = 0.14 ± 0.05 Perryman ct al. 1998; [fe/H] 
= 0.13 ± 0.01 Paulson, Sneden & Cochran et al. 2003). Is 
this discovery a confirmation of the primordial mechanism ? 
We believe that is could be a premature conclusion because 
there is always the possibility that this star when was at 
its pre-main sequence phase, could have been highly metal- 
lic enriched by accretion and that this metallicity was after 
diluted by convection during the successive sub-giant and 
giant phases. 

Independent of the presence of planets, Gonzales (2006) 
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suggests that studies searching for trends between [Fe/H] 
and spectral types, as those realized by Dotter & Chaboyer 
(2003) on the Hyades stars must be extended to other clus- 
ters. We can add that this extension can also be applied to 
studies searching contaminated stars in a cluster by their 
position in a H-R diagram as those realized in the Hyades 
by Quillen (2002). 



8 CONCLUSIONS 

With the purpose to explore a possible mechanism of stellar 
self- enrichment by an injection of rocky type planetesimals, 
we realized N-body simulations in which a forced internal 
migration of a planet interacts with supposed planetesimals 
distributed in a plane disc. Several possibilities are consid- 
ered in respect to the planets eccentricities, migration rates 
and mass. Contrary to QH00 where they considered specially 
the effect of 3:1 and 4:1 mean motion resonances between the 
planetesimals and the planet, we found that the 2:1 mean 
motion resonance is the main mechanism to drive planetesi- 
mals to the surface of the star. This mechanism is essentially 
efficient for slow migration rates of (r = 10 5 — 10 6 yr) and 
low planet eccentricities. 

QH00 restricted themselves to study resonances distant 
enough from the planet, since planetesimals with orbits that 
become planet-crossing are more likely to be ejected from the 
system rather than impact the star. However, the 2:1 reso- 
nance has a protection mechanism that avoid close encoun- 
ters between the planetesimal and the planet. Then plan- 
etesimals locked in the 2:1 resonance are allowed to reach 
high eccentricities and impact the star. The case of e p = 
was also not considered by QH00 because for capture in the 
4:1 and 3:1 resonances it is necessary a minimum eccentric- 
ity. Nevertheless, our results (Figure 6) show that for e p — 
almost all planetesimals are captured in the 2:1 resonance 
and impact the star along the planets migration 

Considering the necessary discs masses in order to pro- 
voke a migration, we are able to set the most realistic con- 
ditions for an enhance of the metal content of the stars. A 
migration of a one Jupiter mass planet from 5 Au up to 
0.01 AU requires excessively large disc masses. The mass re- 
quirements are 10 times smaller, and became realistic, for a 
migrating 0.1 Jupiter mass planet. Our simulations with this 
last kind of planet and for a slow migration rate (r = 10 6 
yr) stopped at pericentric distances between 0.19 and 1.79 
AU when the last planetesimals are removed from the sys- 
tem. Taking into account the mass of the total impacting 
planetesimals we obtain maximum metal enrichments of the 
order of [Fe/H]=0.18 dex. produced at ages of the stars be- 
tween 20 and 30 Myr. These results need however, primitive 
disc masses three times larger than that of a primitive solar 
disc nebula. These calculations open the possibility never- 
theless, width larger primitive discs (but possible to exist), 
to explain the metal abundances of the observed hot Jupiters 
exoplanets located at 0.05 from the stars. 
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